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Abstract: 

We investigate the high-energy scattering in the spontaneously broken Yang - 
Mills gauge theory in 2+1 space-time dimensions and present the exact solution 
of the leading In s BFKL equation. The solution is constructed in terms of special 
functions using the earlier results of two of us (L.N.L. and L.S.). The analytic 
properties of the t-channel partial wave as functions of the angular momentum 
and momentum transfer have been studied. We find in the angular momentum 
plane: (i) a Regge pole whose trajectory has an intercept larger than 1 and (ii) 
a fixed cut with the rightmost singularity located at j — 1. The massive Yang - 
Mills theory can be considered as a theoretical model for the (non-perturbative) 
Pomeron. We study the main structure and property of the solution including 
the Pomeron trajectory at momentum transfer different from zero. The relation 
to the results of M. Li and C-I. Tan for the massless case is discussed. 
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1 Introduction 



Recent experimental data from HERA JTJ on deep inelastic scattering at small x 
and fixed Q 2 and from the TEVATRON on high-energy diffraction revived the 
interest in the longstanding problem of the Pomeron structure and of the relation 
between soft and hard processes at high energy. 

For the hard Regge processes one can use the BFKL theory ||, but we are 
lacking a selfconsistent theoretical approach to the soft Pomeron and have to 
rely merely on general properties of analyticity, causality and crossing symme- 
try in developing an extended and successful phenomenology of high-energy soft 
interactions J§, §, §. 

Some theoretical understanding of the Pomeron has been derived from the 
study of the leading Ins approximation of superrenormalizable models like A0 3 
in 3 + 1 dimensions. The main features of the result have been included in 
the parton model of peripheral interactions and they are the basis of our under- 
standing of the Pomeron structure 0, @. However, such models result in Regge 
singularities with intercept around —1 and do not reproduce essential features of 
the Pomeron. Much effort has been applied to show the selfconsistency of the 
Pomeron hypothesis in the framework of reggeon field theory or Gribov's Reggeon 
calculus ||. A Reggeon field theory approach to QCD has been developed in jlTJ. 

On the contrary, for the hard Pomeron we can apply perturbative QCD and 
derive a number of detailed predictions |TT]]. A special role plays the BFKL 
pomeron J| appearing in the leading Ins (~ In-) approximation. The main 
features of the BFKL Pomeron, however, look different from properties of the 
soft Pomeron. 

In this paper we study the BFKL Pomeron in spontaneously broken 2 + 1 
dimensional gauge theory, using previous results obtained in ref. |12|]. One can 
consider this theory as a simple model for the soft Pomeron. Indeed we show that 
the resulting BFKL Pomeron is a normal moving Regge pole with its intercept 
op(0) > 1. 

The coupling of this theory has the dimension of mass. The interaction is 
superrenormalizable. This results in the absence of scaling violations of structure 
functions due to ultraviolet divergences. On the other hand the infrared singu- 
larities in the massless limit are stronger compared to 3 + 1 dimensional QCD. 
The comparison allows us to discuss the influence of the ultraviolet and infrared 
singularities on the Pomeron structure. 

In QCD (massless gluons in 3 + 1 dimensions) the known way |13|] of solving 
the BFKL equation relies on conformal symmetry. This approach is useless in 
the case of massive gauge bosons. Up to now the solution is not known for the 
massive case. In the special case of 2 + 1 dimensions, however, the equation 
exhibits a simple iterative structure which allows to construct a solution. The 
experience gained in the 2 + 1 dimensional theory will be helpful in solving the 
corresponding equation in the physical case. 
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We obtain the exact solution both for the forward and non-forward cases, 
and calculate the partial wave amplitude for the scattering of two massive gauge 
bosons. We investigate the Regge singularities in the complex angular momentum 
plane and their behaviour in dependence of the momentum transfer. 

The paper is based on an early investigation by two of the authors [12], where 
the basic idea of the iterative solution was formulated for the general non-forward 
case. This investigation was motivated in particular by ||15| , where the BFKL 
equation with the infrared regularization has been considered. We discuss the 
relation of our result with the one by M. Li and C.-I. Tan jnj where the massless 
2 + 1 dimensional gauge theory has been considered. 



2 BFKL equation with massive gluons 
2.1 3+1 dimensions 

Let us start with the short reminder of the results obtained within the leading 
logarithmic approximation of perturbation theory (LLA) for the amplitudes of 
the high energy scattering in the spontaneously broken Yang-Mills theory 0]. 
We discuss the simplest case of SU(2) gauge group with symmetry breaking by 
one Higgs doublet (fundamental representation). This is the case discussed in ||; 
we shall follow the notation of that paper. The generalization to SU(N) gauge 
group is straightforward and is done in the section 3.3. Notice that the details 
depend on the way of symmetry breaking. We consider the case that all gauge 
bosons become massive. 

The amplitude describing the elastic two-particle scattering AB — > A'B' can 
be decomposed into the amplitudes with definite isotopic spin T in the t— channel, 
with T = 0,1,2. 

AiE' = t aa ,a^t bb , + r AA ,A^r BB , + n A ,A^v% BI (2.1) 

The constants V in ( |2.1| ) depend on the kind of scattering particles (gauge 
bosons, fermions, higgs particle), and they are all proportional to coupling con- 
stant r oc g, for their explicit forms see 0. 

In what follows we shall concentrate ourselves on the singlet part of the am- 
plitude (|2.1|). A(°) is related to the partial wave F u) (q 2 ) in the following way 
(j = l+u): 

(5+ioo _ ■ 

A W (s,q) = ^ f du(^) ( —^F^), t = -q\ (2.2) 

<5— ioo 

Kiq2) = Ijty I W^W=W^) UK q " k)Mq2) ' (2,3) 
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A (q 2 ) = -2(g 2 + \m 2 ) . (2.4) 

We write here and in the following the scalar products of transverse momenta 
in Euclidean notation. The function f w {k,q — k) satisfies the following integral 
equation (see Fig.l for notations and graphic form of the equation). 



[uJ-a(k 2 )-a((k-q) 2 )]Uk,q-k) 



to 



A^'S? 1 ^ ^)((1T- q) 2 + m 2 ) K ^ **' q)UK q ~ h) (2 ' 5) 
with the kernel 

K(k,k u q) = A (q 2 ) (2.6) 
+ 7T ttt— — 1 \{k 2 + m 2 )((A;i - qf + m 2 ) + {k 2 + m 2 ){{k - q) 2 + m 2 )] 

and the Regge trajectory of the masive gluons 

/.2n ^ 2 (A: 2 + m 2 ) /• rf 2 ^ 

1 J J (2tt) 3 J (k 2 + m 2 )((fc 1 - k) 2 + m?) ' 1 ' ' 



As a consequence of the integral equation ( |2.5| ) it is possible to express the 
partial wave F^q 2 ) through the solution of Eq. Q2.5|) on the mass-shell 

fw{k, k - q)\ k 2 =iq _ k) 2 = _ m 2 = F w (q 2 ) + A^\q 2 ) • (2-8) 



2.2 2+1 dimensions 

In Refs. [jl2[ p~5f| it was established that in 2+1 dimensional space-time the high 
energy scattering amplitudes derived in LLA are given by formulae similar to the 
ones from the previous section. The main difference is that the transverse space 
in this case is one-dimensional, so the substitution 

d 2 k dk 
(2vr) 2 2^ 

should be made. 
Therefore 

F ^ q2) = uj^f I (k 2 + m 2 )((k- q y + m 2 ) Uk > q ~ ' (2 " 9) 

where now g 2 carries the dimension of mass. 
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Figure 1: The graphic form of the BFKL equation. 



The Regge trajectory for massive vector bosons in the 2+1 dimensions is given 
by the simple rational expression 



g 2 (k 2 + m 2 ) 



a(k 2 ) = 



(2n)< 



g 2 k 2 + m 2 
{k\ + m 2 ){ki — k) 2 + m 2 ) lixm k 2 + 4m 2 



. (2.10) 



For the function f u (k, q — k) we have here the one dimensional Bethe-Salpeter 
type equation 



[u-a(k 2 )-a((k-q) 2 )}Uk,q-k) 



+ 



jrK(k, hi, q)fu,(ki, q - h) (2.11) 



A(q 2 ) (2ir) 2J (ki + m^dh-qy + m 

Our aim is to find the analytic solution of this equation. It is convenient to 
introduce the dimensionless variables 



a s = 



9_ 

4:11171 ' 



e = 



9* 
2nmuj 



2as 



, k — > mk , q — > mq . 



(2.12) 
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Then Eq.( p.ll[ ) takes the form 



fc 2 +4 



+ 



(fc-g) 2 +l - 
(fe-g) 2 +4' 

2 



fa>{k,q- k) 



A' 1 



,(fcf+l)((fci-g) 2 +l) 1 (fc-fci) 2 +l 
An 



fc 2 + l I (fc-q) 2 + l 



-(2g 2 + 



(fci-g) 2 +l 
2/ 



J Z7T 

UkuQ-h) (2.13) 
(2.14) 



The remaining equations of the previous section are unchanged. 

In order to present the main steps of our method for finding the exact solution 
of Eq. (|2.13|) we consider first the simpler case with vanishing momentum transfer 
q = 0. 



3 Forward scattering at high energy. 

In the case with vanishing momentum transfer q = 0, the Eq. fl2.13| ) takes the 
simpler form 

(l + 2e)(A; 2 + A^ 



k 2 + A 



~fu>{k) — A 1 



+ 77 — UvT—-, ,.2 , \ f^ k v > I 3 - 1 ) 



2tt V(^ + i) 2 (k-hy + lkf + l 

where we have introduced the convenient notation 

A* = *±* . (3.2) 
1 + 2e V ; 

In this case we present the methods of solution both in coordinate and in mo- 
mentum representation. In this way different aspects of the problem will be 
illuminated. 

3.1 Coordinate space analysis 

We find that it is convenient to work with the function ^(k) = (k 2 + l) _1 / a; (/c). 
Equation ( |3.1| ) is a linear inhomogeneous integral equation. We try to solve it in 
coordinate space by introducing 

(f) u (k) = J dxe ikx ^{x) . (3.3) 

The main advantage of the coordinate space is the fact that the BFKL kernel in 
Eq. flOp looks simple due to the relation 

*j£_=i e -H. (3.4) 
2tt k 2 + 1 2 y ' 
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Substituting Eq. ( |3.3| ) and Eq. ( |3.4| ) in Eq. ( |3.1| ) we obtain 



(1 + 2s) ( + A 2 ) = ^-(3 e~ N + 25(x)) (3.5) 

a z x 2A 

+2e(~ + 4)e-l^(x) - ^(3e~^ + 2%)) / dy<j>M e~M 
d z x 4 J 

where 5(x) is Euler 5 - function. We shall analyze Eq. fl3.5Q without the inhomo- 
geneous term in order to investigate the leading eigenvalue e. 

<f> u (x) should be bound for the Fourier transform (|3.3j ) to exist. At large 
\x\ only the left-hand side of (|3.5|) is important which leads to the asymptotic 
solution e~ x ' x '. 

Clearly the solution depends on |x| only because the kernel K{x) is an even 
function of x. We introduce a new function 

* u (x) = [1 + 2e(l - e-W)]0 w (z) . (3.6) 

and a new variable z = e _ ' x '. 

For the function <& w the equation looks as follows 

d d$ u (z) „ d$Jz) 13 . , , /0 _s 

dz dz dz 

®Jz) A n e _ r 1 . . §Jz') 

6e 



+ ^(3z + 2S(z-l)) f 1 dz' - 
4 Jo 1 



l + 2e(l-z) 4 v v "h l + 2e(l-2') 

Comparing the coefficients in front of 5- functions we obtain 

d<$>Jz) , eAi f 1 , , 



cfc 12-1 4 Jo 1 + 2e(l - z' 

which will give the equation for the position of the pole in angular momentum 
plane (the intercept of the Pomeron) as will be shown below. The second condi- 
tion 

$ w (z) -> at z -> (3.9) 

follows from the large |x| behaviour of w (x) discussed above. 

The important observation is that solution of Eq. ( |3.7Q obeying ( |3 .9|) can be 
found in the form 

$ u {z) = Cz[l + 2e(l-z)} + ^(z), (3.10) 
where $^ 9 (-z) is the solution of the homogeneous equation 

+ = fe i + £i-,) ' (3 - u) 



^ — V > ( 3 - 8 ) 
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and C is a constant in z, which however depends on the function 

C = / dz' - r . (3.12) 

8 + \e Jo 1 + 2e(l-z') K J 



The solution of the homogeneous equation (|3.11| ) can be easily found. We 
obtain, 

<S> h J(z) = Nz x 2 F l (2 + X,-2 + X,l + 2X,^-) = (3.13) 

(1 + 2efl - z))z x . , , 2ez . 

N " 2 F 1 (3 + A,-1 + A,1 + 2A, 



l + 2e ^ ±v ' ' ' l + 2e' 

Here 2-^1 is the Gauss hypergeometric function and the constant iV can be defined 
from the normalization. 

To find the value of e which corresponds to the bound state we have to solve 
Eq. ( |3.8j ) which using well known properties of the hypergeometric function can 
be reduced to the form 

c _ 9 \ 2 4 _ \ 2 

T -^ 2 F,(3 + A,-l + A,l + 2A,^— )+ (3.14) 

(^Zi) 2fl( 3 + A,A, 1 + 2A,1^) = 

2 1 4 — A 2 

' 2 F 1 (2 + A,-1 + A,1 + 2A, 



1 + f £5 2 + A— ' '•' - ' ■■' ' 3 

+ T^A 2Fl ( 1 + A '- 1 + A ' 1 + 2A '^-)) 

We solved this equation numerically and obtained the value e = e = 4.5934 
which leads to the rightmost singularity at u = u = = 0.436 as in 

accorance with ref. ITS] . 



The way we have solved the BFKL equation is reminiscent of the standard 
procedure of calculating bound states. The rightmost singularity in u, a pole, 
corresponds to the ground state. In the following subsection we solve Eq. (|3.1| ) 
in momentum representation. 

3.2 Momentum space analysis. 

We have to solve the linear inhomogeneous integral equation ( p.l|) . It is possible 
to construct the solution directly by iterations. We rely on the relation 

OO 

dk' 1 1 A + l 

2tt (k - k') 2 + lk' 2 + A 2 2A[A; 2 + (A + l) 2 ] 1 ' ' 
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This means that the action of the kernel on (k 2 + A 2 ) -1 can be expressed by the 
shift of the pole position A — > A + 1. 

Let us formally consider the right-hand side of Eq. (|3.1|) as a perturbation. 
We will consider first the solution of this equation in the interval e G [—1/2, oo], 
where A 2 > 0. We obtain the solution for A > first and continue then analyti- 
cally to the complete complex plane in e or u. 

If we omit the right-hand side of Eq.( p.l|) (the zeroth iteration), the solution 

is 

f^(k) can be represented as the sum of the constant term plus the pole 
term ~ ^ 2 jj^2 • In order to find the contribution arising from the next iteration 
(/ w = f(°) + /« + . . .) let us substitute (gig) into the right-hand side of Q. 
We use now Eq. ( |3.15| ) and obtain 

dk\ 



A (\ + 2) 2(k 2 + l) 
4A(A + 1) 2 A(fc 2 + 4) 



I A + 3 



A + l k 2 + (X + l) 



(3.17) 



We write the resulting first iteration /W as a sum of pole terms in k 2 . In this 
expansion there are three terms: the constant term, the pole term ~ ^,2 j|_^2 and 
as a new term, not encountered in the zeroth iteration f£\ the pole term ~ 
fc2 + (j v+1 - ) 2 • The same procedure can be applied to the subsequent iterations. It is 
easy to see that the expansion for the n-th iteration will be given by the sum of 
the constant term plus the pole terms k ^ x i , A^ = A + k — 1, k = 1, . . . , n + 1. 



Therefore it is natural to look for the solution of Eq. ( |3.1| ) in the form [TJ 



OO A 

^) = /» + El^ - A n = A + n-l. (3.18) 

n=l fc + A n 



Let us substitute this ansatz in the equation ( |3.1|) . Comparing coefficients of 
the pole terms on both sides we find the condition 

A n _ 2e (A + n)(A + n + l) 



A n _! l + 2e(n-l)(2A + n + l) 
This recurrence relation has the following solution 

A -A ( 2£ V' 1 (Ai + 2) B - 1 (A 1 + 3) w _x 

An ~ Al [TTYe) (n-l)!(2A 1 + ' (3 ' 20) 
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where (a) n = a(a + l)(a + 2) . . . (a + n — 1). In this way we arrive at generalized 
hypergeometric functions p+iF p Q^"^ +1 \y) [[U|. In particular, the ansatz (|3.18| ) 
leads to 



/«(*) = fo + ijr^^&iS£SiA\v) ■ (3-21) 



k 2 + X 2 ' 

with A n = A + n — 1, A = Ai and y = Y^Ze' 

There are still two coefficients fo and A\ undetermined in our solution (|3.21| ) . 
The information contained in equation (|3.1|) which has not been used yet can be 
expressed in terms of two conditions. The first condition appears as a result of 
the comparison of residua of the pole term ~ ^2_|_^2 ( the pole at k 2 — > X 2 has 
to be considered separately from other pole terms ~ fe2 \ 2 , n ^ 1). The second 
condition appears as a result of the comparison of the constant terms appearing 
in the expansion, or in other words, considering the left and the right-hand sides 
of Eq. Q at k -> oo. 

Therefore the coeficients fo,Ai are the solution of the following inhomoge- 
neous system of linear equations 



— fo ' a ll + Ai • Cli2 



^7 = fo ■ a n + A 1 ■ a 22 , (3.22) 



with 



A ^ 4(A 2 - 1) 
an = — + 



4 (A 2 - 4) 



„ _ _ 6 , Ag\ 3 771 /A4,A4,A,A2 I 

"12 - ( A 2- 4 )2 "T 4AAj 4r 3l2A+l,A3,A 3 l 
4 rp /A4,Aa,A— li \ , 2A2 771 /A4,A3,A— 1,1 



^a^ffiri!/) + ^^&l:tr\y) (3-23) 



_ A 2(2A 2 - 5) 
° 21 " T + (A 2 - 4) 



_ A) A3 p /A 4 ,A 4 ,A,A 2 I \ I 2 , A4iA . 

^22 - |^24^3V2A+1,A3,A 3 I2/J + l\2\+l\ 



Using well known relations among the hypergeometric functions [16] it is 
possible to express all higher hypergeometric functions through the two basic 2 -Pi 
functions 

fa = 2^1(2^1 1 v) 



h = aftG&ilv) (3.24) 
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We quote here only one of these relations 

Wittily) = f« { x-i )x % + fb (\-iyxixi ■ (3 - 25) 

The solution of the system ( |3.1| ) expressed in terms of the functions f a and 
has the form 



fc 



AjA§(A-2) 

(34A 2 -64)(A 2 
1 A(A-2)A 3 



4 



(fa 



(13A 2 -16) 
A 



-J) , f 48A 2 -84 
•/ 6 (A-2)A 3 



+ 18/ 6 ) + 
i 



3A t 



(A-1)(A-2)A§A: 



:(/„A(l + 2A 2 ) + / ( 



3(A 4 + 2) 
(A - 1)A 2 



(3.26) 



These formulae together with Eq. ( |3.21| ) represent the solution of the integral 
equation fl3.1|) . 



It should be noted that in Ref. [12] instead of the first equation of the system 



fc 2 +A 2 



( |3.22|) (resulting from the comparision of the residua of the pole terms 
appearing on both sides of eq. ( p.l| )) another boundary condition was used, the 
absence of the normal thresholds, 



Uk 2 ^-A)=2Uk 2 



-I) 



(3.27) 



This condition can be derived from the equation ( |3.1| ) if one requires that fu,(k) 
is a regular function in the neigbourhood of k 2 = —4. In terms of our ansatz this 
condition reads 



fo = A 1 



1 



A 2 



'A 3 ,A 3 ,A-2 
v2A+l,A 



l\y)- 



A 2 



"3^2 



\2,A-1 

>2A+1,A 



\y) 



(3.28) 



The iterative solution of the equation (|3.1|), / w (/c), as described above, is a 
function which is by construction regular in the points k 2 = —n 2 . Therefore, the 
condition ( |3.28 ) should not give an additional restriction on the function (3" 



as compared with the conditions given by the system ( p.22| ). Indeed, expressing 
the hypergeometric functions in (|3.28| ) in terms of the functions f a and it can 
be checked directly that the difference of the two equations in ( |3.22| ) and the 
condition (|3.28| ) are equivalent. 



3.3 Regge singularities of the forward partial wave 

We discuss now the implications of the obtained solution fuj(k) for the partial 
wave of the scattering amplitude F u . The partial wave can be calculated 
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either by ( |2.9| ) or by the mass-shell relation ( |2.8| ). We have checked that the both 
methods lead to the same result 

F «(9 = ) = TTT7' (3 ' 29) 

where 

/ = _1_ ( (34A--64)(^-l) /a + A(48 , 2 _ g4) j ( (13^-16) A + ^ 

(3.30) 

Let us discuss the singularities of F w considered as a function of the complex 
variable uj. The hypergeometric functions are defined in terms of the hypergeo- 
metric series which are convergent inside the circle of unit radius in the variable 
y = 2e/(l + 2e). The continued hypergeometric functions are analytic in the 
complex plane of their argument y, with a cut from y — 1 to y — oo. In the e 
plane this corresponds to the cut appearing on the interval e G [— oo, — 1/2]. 

As a function of their parameters a p+ i, (3\,...(3 P the hypergeometric 

functions have only simple poles if one of the lower parameters fii, . . . f3 p ap- 
proaches a non-positive integer value n. We see from (|3.2|), ( |3.24j ) that both f a 



and ft have poles of this origin at A = yf^ff — > ~^ L - Therefore these poles lie 
on the second (unphysical) sheet of the square root. 

Further singularities of and, consequently, of F^ appear at points, where 
the determinant of the system of linear equations ( |3.22j ) vanishes, i.e. at the zeros 
of the denominator in (|3.29| ) 

A + f = 0. (3.31) 

This results in poles in uj. 

Analyzing the condition (|3.31| ) numerically outside the interval e G [— oo, —1/2], 



where the cut is located we have checked that there is only one Regge pole in the 
vicinity of the real axis located in 

Q 2 

e = e = 4.5934 , u = u = — . (3.32) 

27rmeo 

The result coincides of course with the one obtained in the coordinate represen- 
tation. Therefore we can conclude that at q = the partial wave F w has the 
following singularities on the physical sheet of complex u plane. There is a finite 
cut on the negative part of the real axies covering the interval uj G [to , 2 ,co , i], with 

2 

to>2 = iO\ = 0. And there is a single pole in the positive part of the real axis 
at uj = u>o, ( |3.32| ). 

Let us discuss the nature of the singularities at the branch points. Near the 
right end-point of the cut, uj = uj\ = we have e — > +oo, A — > 1, y — > 1 and 

/„ = -2(1 + log ^^) + 0((A - 1) log(A - 1)) , 
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/ 6 = 2 + 0((A-l)log(A-l)) 
Therefore the partial wave behaves like 



(3.33) 



F u -> -A~ l + IQAf/ log(A - 1) (3.34) 

Near the left end-point of the cut, to = uj 2 we have e — > —1/2, A — > +oo, 
y ~ — A 2 /3 — > — oo and 

A(logl2-21ogA) 

/a = = (1 + 0(1/A)) (3.35) 

VttA 

A(logl2-21ogA) A1 A 

A = " ^ 6 (1 + 0(1/A)) (3.36) 

Therefore the partial wave behaves like 

1 1 

FuJ "Aq + 24/3 ~ 72- (A + 24/3) 2 log A (3 ' 3?) 

Note that the solution of the corresponding homogeneous equation can be 
obtained from the solution of inhomogeneous equation which we have just 
found. The spectrum consists of one discret level uj = uj and the continuous 
part uj G [a; 2 ,o;i]. The corresponding eigenfunctions can be found as follows: the 
residue of F u of the pole at uj = uj gives (up to the normalization constant) 
the wave function of the discret level and by calculating the discontinuity of F w 
on the cut it is possible to find the eigenfunctions belonging to the continuous 
spectrum. 

We would like to add a comment on how the results depend on the number 
of colours N. In the case of arbitrary N we have to substitute in the above 
equations H 

2 2 N A «, 2 , N 2 + 1 2^ 

9 ^9 Y> A)->~2(<? + N2 m ) 
N 

e-ey. (3.38) 
As in the case N = 2 there is a leading Regge pole at arbitrary iV located at 



UJ, 



(N) 
) 



uj {N) - 9 - - 1 f3 39^ 

^ - 2^ y ^) • (3 - 39) 

is calculated in analogy to e above. We find that decreases slowly 
with N approaching a limit ef^: = 4,5934, e[, 3) = 3.8000, = 3.5693, 
eg* = 3.3025. 
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4 Non-forward scattering 



4.1 Solution of the equation 

The main steps which have been made in the section 3.2 to derive the solution of 
the forward equation can be generalized to find the solution of the non-forward 
equation ( |2.13| ). The expression appearing on the left-hand side of Eq.( |2.13| ) in 
the square brackets [. . .] can be rewritten in the form 

(l + 2e)[x 2 + (A+) 2 ][x 2 + (A-) 2 ] 
L-"J [( x _ g /2)2 + 4 ][(x + g/2) 2 + 4] 1 ' ] 



where 



x = k — q/2 



(A±) 2 = i±^-^± 



9e 2 -g 2 (4 + 5e)(l + 26) ^ 



l + 2e 4 \| (l + 2e) 2 

Now, in analogy with the iterative way of finding the solution for q = , we 
see, that the zeroth iteration can be expanded into a sum of the constant term, 
and two pole terms: ~ 3 ,2+(a+) 2 an< ^ ~ aP+DFj 2 • ^ should be noted that the zeroth 
iteration for f u (k, q — k) depends on the specific combination of the momenta k 
and q — k, i.e. it is a function of the variable x 2 = (k — q/2) 2 . The notation x 
should not be confused with the position. Calculating the next iterations it can 
be seen that this feature remains true and the solution can be represented in the 



following form [12 



At 



A: 



\2 ' n 



A± = A ± + n - 1 . 



(4.3) 



Substituting this ansatz into the equation ( |2.13p we find two recurrence rela- 
tions similar to the one for the case q = 0. Their solutions can be written in the 
form 



with 



\£ _ ( 2e \ 
if Vl+2e; 



n-1 



(n-1) 



X 



(A±) n _!(A± + f ) n _ 1 (A±-f ) n _ 1 (A ± -d+) n -i(A ± -dt) n _i(A ± -d;) n -i(A ± -rf:) n - 1 
(A±)„_i(2A±+l) n _i(A±+AT+l) n _i(A±-A=F+l) n _i(A± + f ) n _i(A±-f )„_i 



(4.4) 



(4.5) 



It is possible to rewrite our ansatz Eq. ( |4.3| ) in terms of the generalized hyper- 
geometric functions 
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Ux 2 ) =fo + 
At- 1 



9-^81 



x 2 + (A+) 2 



f \+,\++^,\+- i ^,\+-d\,X+-dt,\ + -d-,\+-dZ,\ + +ix,\+- 
A+,2A++1,A++A-+1,A+-A-+1,A+ + H a+-H i A++m:,A+-m 



"\v) + 



(4.6) 



The conditions which determine the coefficients fo, Af, A[ are also analogous 
to the ones used in the case of q — 0. 

The condition obtained by taking the limit x — > oo, or k — > oo, has the form 



l 



/o 



An 



+ 



2e-l 



+ 



.+ [ A A+ „ .A+.A+H-f ,A+-f ,A+-d+ A+-dl,A+-d;,A+ 



A+(q2+4)((A+)2 + V) 



2A^ 



A+,2A++l,A++A-+l,A+-A-+l,A+ + ^,A+-f 

" \y) 



A+((A+)2 + ^) 7 6l A+,2A++l,A++A-+l,A+-A-+l,A+ + f,A^ 



)+ 



(4.7) 



We have found that it is convenient to use as the last two conditions the 
absence of normal thresholds (see discussion at the end of of section 3.2) 



2/ w ((|±0' 



(4.8) 



The condition corresponding to the lower signs is 



h = At 



((A+)2-(2+i2)2) 



8-^71 



/ A 2 > A 3 +T' A 4 ' 



,A+-2- 



,A+ -dt,A+ -d+ ,A+ -CA+ -<T 



A+,2A++l,A++A-+l,A+-A-+l,A++f ,A+-f ,A+-l-f 



7-i ' 4 T 2 ' 4 



,A+-1- 



,A+-d+,A+-d+,A+-d~,A+-d- 



((A+) 2 -(l+f ) 2 ) 8 7 ^A+,2A++l,A++A-+l,A+-A-+l,A+ + f ,A+-f ,A+-f 

+^r[A+ <-> a-] . 



|y) 



|y)- 



(4.9) 



The other equation is obtained from the above one by the substitution q <-» —q. 
The difference of these two condition can be written in the limit q — > as 

+ 0(1) j + A\ {-iCq + 0(g 2 )} = (4.10) 



Therefore 



Ai = A^[Cq 2 + 0( q 3 )] , 



(4.11) 



where C is some constant. 
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We see that at q — > the series of poles ~ a , 2+ (^+)2 decouples from the solution 
in accordance with our previous considerations for q = 0. 

In this way we have solved the equation ( |2.11| ) for arbitrary momentum trans- 
fer q. The solution / w is given by (|4.6f) with the coefficients fo, Af and A{ 
determined from linear system of equations (4.7|), (|4.8| - C9 ). 



4.2 Properties of the partial wave 

We investigate the partial wave F^q 2 ) obtained from the solution by Eq. Q2.9|) 



J ° ^ A +[(A+)2 + ^] 7 6l A+ A+ + f ,A+-f ,2A++1,A++A-+1,A+-A-+1 ^ 



6A 

.12 \ — _i2 



A~A^ ,AJ,AJ + H iA J-H,A--d+,A--d+,A-- f Z;,A-- C r , . 

+ A-[(A-)2 + ^] 7 6l A3,A 3 - + H,A 3 --H, 2 A-+l,A++A-+l,A--A++l ^ 



(4.12) 



First of all it should be noted that all equations above are written under the 
assumption that we choose the convention for the square root expression for \^ 
with the real parts of A ± being positive for small q and real e > — |. 

If q 2 > A ± are complex conjugate numbers if e belongs to the interval 

[ei,oo], where €\ = 139 2 (9-wq j y~~ - Since —1/2 < A x < 0, for any positive uj and 
therefore for any positive e, A ± are complex conjugate. Let us choose by the 
definition as A + the expression which has negative imaginary part (for q > 0). 

In the following we study the Regge singularities and the behaviour at — q 2 = 
t — > — oo and at positive t up to the first threshold t = 4. 

Since the solution behaves smooth at q — > we conclude that at small q the 
structure of the Regge singularities is similar to what we have found for q = 0. 
The position of the leading Regge pole depends on t = —q 2 . The result of the 
numerical calculations is plotted in Fig. 2 for values of t from —4 up to the vicinity 
of the first threshold at t — 4. The trajectory is almost linear in the vicinity of 
t = with the approximate slope 0.34^f as shown in Fig. 3 . 

We would like to mention that the Pomeron trajectory has about the same 
slope ( a' P (0) ) as the gluon trajectory. More interesting would be, within the 
same approach, to compare the Pomeron trajectory with the Reggeon trajectory. 
In order to do so one has to calculate the Reggeon trajectory in 2 + 1 QCD using 



the techniques developed in Ref. jn[]. It will be a challenging problem for the 
future. 

At larger \t\ the trajectory deviates strongly from the linear behaviour. It goes 
to infinity for t approaching the threshold t = 4 and returns from — oo above the 
threshold. The behaviour of the Pomeron trajectory near t — 4 has been obtained 
also by solving Eq.(2.13) in the asymptotics of large u and t — > 4 with the result 
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Figure 2: The trajectory of the pomeron pole in units of The momentum 

transfer is given in units of m 2 . 
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t 



Figure 3: The behaviour of the pomeron trajectory in the vicinity of t = 0. 
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This confirms the numerical result of Fig. 2. 

2 

The branch points u)\ = 0, cu 2 = — - 2 — do not depend on t. However the 
singularities located at the unphysical sheet can come up to the physical sheet 
as t increases. Indeed both the poles arising from the lower coefficients fa in the 
hypergeometric functions and the poles at the vanishing determinant depend on 
q. There are also branch points arising from the square roots in the expressions 
of A ± in terms of to ( |4.2|) the position of which depend on q. The numerical 
investigation of the solution shows that besides of the pole, which was originally 
the leading one, another pole emerges from the unphysical sheet if t crosses the 
threshold value. 

Now we investigate the behaviour at t — > — oo. In the limit of large q we have 




and 

d b a = a^ + T ™ , rW = a ^ ~ 1 . 
2 ' 2 

Inserting these relations into the linear system we have found that asymptot- 
ically in q 



h 




where 



k+2 k-2 



f T? /A:+r(+),fc+r(-),fc+2i \ 

Tl - 3^2bifc+l,fc+l W) 

H - 3^2bfc+l,fc+2 W) 

J3 — 4-T3i,2fe+l,fc+l,fc-l \y) 

u — 4-T3(,2jfc+i,fe+i,fe \y) 



As the result we obtain from Eq.(|4.12|) the asymptotics of the partial wave 




k+V y k + l k-l k + 2 k-2' 
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The behaviour of F w near the right branch point oo = is i^l^-nx; ~ 
This is to be compared with the behaviour at the same point for q = 0, F w \ q 2 =0 ~ 
const . The numerical calculation of the Pomeron trajectory ( see Fig. 2 ) shows 
that the pole is moving towards the right branch point with decreasing t. From 
both observations we conclude that the Pomeron pole moving with t reaches the 
right branch point oo = asymptotically for t = —q 2 — > — oo. 



5 Comparison with the massless case 

M. Li and C.I. Tan [|E4j] investigated 2 + 1 dimensional QCD without symmetry 
breaking, i.e. for massless gluons, and obtained just a fixed cut starting at j — 1 
as the leading singularity in the vacuum exchange channel. We try to understand 
the relation of their result to ours, in particular whether the pomeron pole is 
absent in the massless case and how it disappears at m —>■ 0. 

The infrared singularities in 2 + 1 dimensions are stronger compared to 3 
+ 1 dimensions. The limit m — > has to be performed carefully. Clearly, the 
scattering amplitude of vector bosons has no finite limit at m — > 0. We consider 
the scattering of two colour dipoles of transverse sizes xi, X2, which is the case 
studied in ||14| . The partial wave of the dipole-dipole forward scattering is given 



by the convolution of two dipole impact factors |14| (here xq is the size of the 
dipole ) 

$d(x, k) = Asin 2 kx (5.1) 
with the Reggeon Green function 

The Reggeon Green function is the particular solution of the BFKL equation with 
5-functions as inhomogeneous term. It is related to our solution fu(k) which is 
more closely related to the vector boson scattering as follows 

A °° d k 

u(k4m*) Uk) = I G ^ ki) wrb ■ (5 - 3) 

— oo 

Near the pomeron pole we have 

G^y.^MWM, Wo = ^ , (5.4) 

where ipo{k) is the wave function of the two-boson bound state corresponding to 
the pomeron. It is normalized to 1 and can be obtained from f w by studying 
(|5.3|) near oo . e is the number quoted in Eq.(3.32). 
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Restoring the mass dependence we obtain from the solution ( |3.21| ) 



Uk) = — 2 4>(-^-,-) (0.0) 
m z m Zttuj m 



The solution depends smoothly on k and the integral with a bounded function 
$d(X k) exists. Therefore also ipo(k) has these features: 

Mk) = , (5-6) 

<m m 



with a being some numerical constant. The contribution of the Pomeron pole to 
the scattering of two colourless dipoles with sizes x\q and £20 is given at small m 
by the partial wave 

hg mx 2 l0 x 2 20 . (5.7) 



UJ 



— 3l 



Here b is some number. This leading contribution to the forward scattering of 
of dipoles does not behave smoothly at m -* 0. The pole goes to plus infinity, 
resulting in a divergent contribution. Expanding in g 2 we observe that the diver- 
gence starts at the g A term, corresponding to s-channel intermediate state with 
two additional gluons. 

This observation is confirmed by calculating G^(ki, k 2 ) iteratively and evalu- 
ating the corresponding contribution to the dipole scattering partial wave Eq.(5.2) 
in the following way. We have to iterate Eq.(2.13) with the inhomogeneous term 
replaced by 5(ki — k 2 ), which is the zeroth approximation of G^. Unlike above in 
sect. 3.2 and 4.1 the iteration now proceeds order by order in g 2 or e. Replacing 
Guj in (5.2) by 5{k\ — /c 2 ) we obtain that the region of k%, ki ~ m gives a negligible 
contribution for m —>■ 0. Taking the first order approximation in e for G^ leads 
to a finite contribution of that small- k region. With the 0(e 2 ) approximation 
for G^ we obtain a like — diverging contribution. Starting from this order of 
perturbative expansion the amplitude of forward dipole - dipole scattering does 
not exist in the massless limit. 

Consider now the scattering at non-vanishing momentum transfer. Let us fix 
the value t p h ys in physical units (GeV 2 ) and look at the relation to our dimen- 
sionless variable t = —q 2 (in units of m 2 ) 

tphys =tm 2 . (5.8) 

Provided t v h ys < 0, the corresponding value of t approaches — 00 at m — > 0. Thus 
the pomeron pole approaches the branch point at j — 1. 

The singular contribution ( |5.7| ) appearing only at t phys = is absent in the 



infrared finite dipole scattering amplitude constructed in [14 
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Let us now study the massless limit directly in the equation. We restore the 
masses in Eq. Q2.13| ) and do the shift k — > k — | as in Eq. ( }4.1|) 



1\2 +m 2 



(k - |) 2 + m 2 

(£;-§) 2 + 4m 2 ' (k + §) 2 + 4m 2 ' 



fu(k,q) 



A l + em J 



dki 
27 



A 



((* 1 + |)* + m ')((* 1 -|) a + m a ) 



+ 



(A; — ki) 2 + m 2 



(fc + f) 2 + m 2 (A;-f) 2 + m 



+ 



(h + iy + m 2 (k 



Q\2 



+ m 2 



U(h,q) . (5.9) 



We perform the Fourier transformation with respect to k 

dk 



-ikx 



2vr 



fu(k,q) 



(5.10) 



and obtain 



(l + 2e - 2ee- m '*') / w (x, g) - ^em | rfy / w (y, g) cos |(y - x) e -*»l— wl 



— S(x) + eS{x) 
2A em 8(x) 



An 



q 2 + 4m 2 



+ 2 



q(q 2 + 4m 2 ) 



—em e 



dyf UJ (y,q)e mlyl sm^\y\ 



dyfu>(y,q) cos|(y-:c)e m|a; 5/1 (2sgn(x) sgn(x - y) + 1) . 



^^(2/,g)e- m ^ cos|y 
9 1 



(5.11) 



It should be stressed that Eq. (|5.11|) is the general BFKL equation for 2 + 1 
QCD in the coordinate space at any value of the momentum transfer t = —q 2 - 
As discussed above the behaviour at m = is different for forward and non- 

+ 2 



q 2 +4m 2 



forward cases. Indeed the coefficient of second term on r.h.s. e 

vanishes at m — ► for q ^ but behaves like — if we put q = before taking the 
limit m — > 0. We discuss in the following the massless limit in the non-forward 
case. We approximate Eq. (|5.11|) at m — > expanding in particular e - " 1 ^' keeping 



terms em (because e 



27T771 I 



In this way we obtain 



(1 + 2em|x|) f u (x,q) - -em J dyf u {y,q) cos-(x - y) 

-t-6(x) -4=8{x) — f dyfu{y,q) sin^|y| 
/±o q J a 

-2emsgn(x) J dyf u (y, q) cos |(a; - y) sgn(x - y) . 



In terms of the function fu(x, q) defined as 

' d 2 ~ 2 
dx 



fu(x,q) 



(5.12) 



(5.13) 
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Eq, ( |5.12| ) has the following simple form 




1 + 2em\x\)f u (x,q) 



5{x) 



(5.14) 



The solution has the form 



1 + 2em\x\)f UJ (x,q) = 

1 • ?| I , a • 9 , D 1 
= — sm - \x \ + A sm -x + B cos -x 

qA 2 1 1 2 2 



(5.15) 



where A, B are some arbitrary constants. 

The original function fuj(x,q) can be calculated readily. We write here only 
part of the result with A = B = 0, which corresponds to certain boundary 
conditions 



1 



5{x) 



2em cos %\x\ 



+ 



(em) 2 sin 



l + 2em|x| (l + 2em|x|) 2 g(l + 2em\x 



(5.16) 



This result has similarities to the expression for the dipole density n w (xo,x,g) 
derived by Li and Tan (Eq.(3.4) in the second paper of Ref. fl4fl ). Our fuj(x,q) 
is not the dipole density and our equation does not know about the dipole size 
on which n^^xo, x, q) depends essentially. However, introducing the dipole size 
Xq > by replacing the r.h.s. of ( |5.14|) by — 5(xq — \x\) and restricting the range 
in x to \x\ 7^ 0, we obtain 



G u (x,XQ]q) = 

5(xo — \x\) 
1 + 2 em \x\ 
4em cos | (xq 



d " L ^ 1 2 at I u • 9/ 
^ + T } [l + 2em|x|] / (3:o " N)sm 2 (Xo 



N) 



+ ^(x - |x|) 

- kl)' 



16(em) 2 sin |(x — |x|) 
g(l + 2em|x|) 3 



;i + 2em|x|) 2 



(5.17) 



We denote by G w (i, x ; (?) the analogon of /L,(x, g) ( |5.13|) of the modified equation. 
This particular solution of the modified Eq. Q5.14D reproduces the dipole density 
n^ixo, x, q) of Ref. [14| up to terms proportional to 5(x). 



6 Summary 

The reduction of the dimensionality to 2 + 1 simplifies the high-energy scattering 
amplitudes and in particular the BFKL equation. The equation can be solved 
analytically even in the case with masses introduced by spontaneous symmetry 
breaking. 



24 



In the forward case we have discussed the solution both in coordinate and in 
momentum space. In the coordinate space the similarity of the Pomeron pole to 
a two gauge boson bound state has been emphasized, whereas in the momentum 
representation the iterative structure becomes transparent, which has been used 
further to solve the equation in the non-forward case. 

We obtain the partial wave for the scattering amplitude of vector bosons in 
an analytic form. This allows us to study the leading and non-leading Regge 
singularities and their dependence on the momentum transfer both for negative 
and positive t. 

At small momentum transfer we find a leading Regge pole, the Pomeron, with 
an intercept above j = 1 and a trajectory approximately linear for small t. Beside 
this pole there is a branch cut whose right end is located at j = 1 independent 
of t. The Pomeron pole approaches the cut for large negative t. 

We have compared our result with the one by M. Li and C.-I. Tan WM and 
have discussed the peculiarities of the massless limit. This limit is different for 
the forward and non-forward case. The massless limit of our result for the non- 
forward amplitude is close to the result found by Li and Tan who used a quite 
different approach. A modified equation (by introducing the dipole size as an 
additional parameter) reproduces the result of these authors. However, in the 
forward case the Pomeron pole leads to a divergent contribution, which is absent 
in the amplitude of Ref. [|14[ . 

Although the model lives in unphysical 2 + 1 dimensions some features can 
be related to the phenomenology of high-energy scattering. The leading pomeron 
pole which is clearly separated from non-leading singularities corresponds to the 
common idea about the soft (non-perturbative) pomeron. Furthermore, a situa- 
tion with a pole with intercept larger than 1 and a fixed cut just at 1 would result 
in a change of the s-dependence with t . It is interesting to note that the constant 
term in the high energy asymptotics, which correponds to the cut contribution, 
has been used to describe the experimental data on high energy behaviour of 
total cross sections |T7|, inclusive spectra [ITS] and diffractive dissociation flI9| . 

The divergence of the trajectory at t — > 4m 2 seems to exhibit an infinite 
series of bound states of two massive gluons. This would differ clearly from the 
features of the hadronic reality. We understand that this is an artifact of the 
leading Ins approximation, since a potential of finite range created by massive 
boson exchange cannot have an infinity of bound states. 

The simplicity of the model makes it useful for further investigations. In- 
cluding fermions the amplitudes with quantum number exchange could be con- 



structed. There will be no direct analogy neither to DGLAP [20| equation nor 



to the nonlinear double-log equation [21]. More interesting could be the study 
of amplitudes with multiple exchange of reggeized gluons, in particular with the 
exchange of negative charge parity (Odderon). 

The model can serve as a testing ground for the non-perturbative treatment 
of diffractive processes [22], [ 23] , and the high parton density effective action [|24{] . 
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Also, the effective action of high-energy scattering J25| and Gribov's reggeon field 
theory ||] can be studied in the simpler situation of 2 + 1 dimensions. 
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